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Abstract 

Thermal  simulation  was  applied  to  2  Wh-class  cells  (diameter  14.2  mm,  height  50  mm)  using  LiNio.7Coo.3O2  or  LiCo02  as  the  positive 
electrode  material,  in  order  to  clarify  the  thermal  behavior  of  the  cells  during  charge  and  discharge.  The  thermal  simulation  results  for  the 
2  Wh-class  cells  showed  a  good  agreement  with  measured  temperature  values.  The  heat  generation  of  a  cell  using  LiNio.7Coo.3O2  was  found 
to  be  much  less  than  that  using  LiCo02  during  discharge.  This  difference  was  considered  to  be  caused  by  the  difference  in  the  change  of 
entropy.  A  250  Wh-class  cell  (diameter  64  mm,  height  296  mm)  was  also  constructed  using  LiNio.7Coo.3O2  and  thermal  simulation  was 
applied.  We  confirmed  that  the  results  of  the  thermal  simulation  agreed  with  measured  values  and  that  this  simulation  model  is  effective  for 
analyzing  the  thermal  behavior  of  large-scale  lithium  secondary  batteries.  ©  2002  Elsevier  Science  B.V.  All  rights  reserved. 

Keywords:  Lithium  secondary  batteries;  Thermal  simulation;  Graphite-coke  hybrid;  Carbon;  LiNio.7Coo.3O2 


1.  Introduction 

Recently,  intensive  efforts  have  been  made  to  decrease  the 
consumption  of  energy  in  Japan,  in  response  to  a  growing 
difference  in  the  electric  power  demand  during  the  daytime  and 
at  night.  To  decrease  the  consumption  of  energy,  one  approach 
is  to  develop  systems  of  efficient  energy  utilization  by  load 
leveling.  The  efforts  of  load  leveling  have  mostly  been  on  the 
supply  side  with  pumped  storage  hydropower  generation, 
and  so  on.  However,  taking  into  account  the  drastic  increases 
in  power  demand  in  urban  areas,  further  efforts  concerning 
load  leveling  are  now  required  on  the  demand  side  by,  e.g.  using 
dispersed  type  battery  energy  storage  systems  for  home  use. 
These  batteries  for  load  leveling  systems  must  be  designed 
to  provide  a  long  cycle  life,  such  as  10  years. 

Lithium  secondary  batteries  using  a  carbon  negative 
electrode  have  a  high  energy  density  and  show  an  especially 
superior  charge/discharge  cycle  performance.  This  is 
because,  in  their  charge  and  discharge  reactions,  lithium 
ions  transfer  between  the  positive  and  negative  electrodes 
and  there  are  no  significant  changes  in  the  crystal  structures 
of  both  electrodes.  Accordingly,  lithium  secondary  batteries 
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using  a  carbon  negative  electrode  are  an  attractive  candidate 
for  dispersed  type  battery  energy  storage  systems.  From 
studies  on  various  carbon  materials  as  the  negative  electrode, 
we  have  found  that  graphite-coke  hybrid  carbon  is  a  pro¬ 
mising  material  as  the  negative  electrode  material  for  battery 
energy  storage  systems  because  of  its  excellent  cycle  per¬ 
formance  and  large  discharge  capacity  [1], 

On  the  other  hand,  in  the  development  of  lithium  sec¬ 
ondary  batteries  for  energy  storage  systems,  the  thermal 
management  of  the  cells  is  very  important.  This  is  because 
these  batteries  are  required  to  have  a  much  larger  capacity 
than  ordinary  consumer  batteries  and  because  their  charge/ 
discharge  currents  become  high.  A  number  of  studies  invol¬ 
ving  thermal  simulation  have  been  made  on  lithium  sec¬ 
ondary  batteries  [2-1 1],  in  order  to  estimate  the  temperature 
changes  during  charge  and  discharge.  Few  studies,  however, 
have  confirmed  a  calculated  temperature  with  a  measured 
value  by  using  actual  large-scale  cells. 

In  this  study,  the  thermal  simulation  of  large-scale  lithium 
secondary  batteries  was  investigated  and  the  simulated 
values  were  confirmed  by  using  250  Wh-class  cells.  The 
thermal  and  electrochemical  characteristics  of  different 
types  of  positive  electrode  materials  (i.e.  LiCoCL  and 
LiNio.7Coo.3O2)  were  also  examined  to  clarify  the  relation¬ 
ship  between  the  thermal  behavior  and  the  entropy  change. 
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Nomenclature 

cp  specific  heat  capacity  (J  g_1  K_1) 

Eq  electromotive  force  (V) 

I,  current  (A) 

Q  heat  generation  in  the  cell  (J) 

r  cell  radius  (mm) 

t  time  (s) 

T  temperature  (K) 

V,  cell  voltage  (V) 

Z  distance  in  the  axis  direction  (mm) 

Greek  letters 

1  thermal  conductivity  (W  m ~  1  K  1 ) 

p  density  (g  cm-3) 


2.  Experimental 

2.1.  2  Wh-class  cells 


2.  The  heat  generation  is  uniform  in  the  cell  during  charge 
and  discharge. 

3.  The  convection  and  radiation  heat  transfer  in  the  cell  are 
neglected. 

4.  The  temperature  dependence  of  the  thermophysical 
properties  excluding  the  specific  heat  capacity  of  the 
separator  is  neglected. 


Under  the  above  assumptions,  the  following  heat  balance 
equation  can  be  written  for  each  part  of  the  equipment. 


dT 

c„P  dt  -  K 


if  T  1  DT 


d2T 


ry  H - 7T~  I  +  ^  ^  Q 


dr2  '  r  dr )  '  "z  dz1 


(1) 


where  cp  is  the  specific  heat  capacity,  p  the  density,  /  the 
thermal  conductivity,  t  the  time,  T  the  temperature,  r  the 
radius,  z  the  distance  in  the  axis  direction,  and  Q  is  the  heat 
generation  in  the  cell. 

All  calculations  were  made  with  explicit  finite  difference 
equations  and  reported  thermophysical  properties  [2,13]  for 
the  materials  in  the  equipment.  The  heat  generation  in  the 
cell  was  calculated  by  the  following  equations  [2,12,13]: 


The  diameter  and  height  of  each  2  Wh-class  cylindrical 
cell  were  14.2  and  50  mm.  LiCoCU  or  LiNio.7Coo.3O2  was 
used  as  the  positive  electrode  material  and  graphite-coke 
hybrid  carbon  was  used  as  the  negative  electrode  material. 
The  graphite-coke  hybrid  carbon  was  prepared  by  mixing 
graphite  and  coke  in  a  4  to  1  weight  ratio.  A  mixed  solvent  of 
ethylene  carbonate  and  diethyl  carbonate  containing  1  M 
LiPFg  was  used  as  the  electrolyte.  A  stainless  steel  vessel 
was  used  for  the  cells.  The  cells  were  charged  with  a 
constant  current  of  65  mA  until  the  cell  voltage  reached 
4. 1  V  and  discharged  with  a  constant  current  of  500  mA  until 
the  cell  voltage  reached  2.7  V. 

2.2.  250  Wh-class  cells 

The  diameter  and  height  of  each  250  Wh-class  cylindrical 
cell  were  64  and  296  mm.  LiNio.7Coo.3O2  was  used  as  the 
positive  electrode  material  and  the  other  materials  were  the 
same  as  those  used  for  the  2  Wh-class  cells.  The  cells  were 
charged  with  a  constant  current  of  8.8  A  until  the  cell  voltage 
reached  4.2  V  and  discharged  with  a  constant  current  of  70  A 
until  the  cell  voltage  reached  2.7  V. 

3.  Model  and  assumptions  of  the  thermal  simulation 


Gtot  —  Gent  +  Gohm  (2) 

Gohm  =  (v,  —  Eq)  X  I,  (3) 

Gent=r(^)x/,  (4) 

where  Gtot  is  the  total  heat  generation,  Gent  the  heat  gen¬ 
eration  due  to  the  change  in  entropy  along  with  the  electro¬ 
chemical  reactions  on  both  electrodes,  Gohm  the  heat 
generation  due  to  the  polarization  including  the  ohmic 
voltage  drop,  V,  the  cell  voltage,  I,  the  current,  and  E0  is 
the  electromotive  force.  In  Eq.  (2),  we  neglect  the  heat 
generation  due  to  the  self-discharge  of  the  cell  and  the  side 
reaction  in  the  cell. 

In  this  study,  the  open  circuit  voltage  was  used  instead  of 
the  electromotive  force.  In  Eq.  (3),  the  internal  cell  resis¬ 
tance  is  treated  as  a  constant  during  charge/discharge.  The 
value  of  ( dEQ/dT)  is  calculated  by  approximate  equations 
obtained  from  the  experimental  results  shown  in  Fig.  5. 

The  positive  electrode  material  was  composed  of  90  wt.% 
LiNi0  7Co0  3O2,  5  wt.%  polyvinylidene  fluoride  (PVdF),  and 
5  wt.%  carbon.  The  negative  electrode  material  was  com¬ 
posed  of  76  wt.%  graphite,  19  wt.%  coke,  and  5  wt.%  PVdF. 

Cell  parameters  and  their  values  for  the  thermal  simula¬ 
tion  of  2  Wh-class  cells  are  listed  in  Table  1 .  The  heat 


For  the  thermal  simulation,  we  followed  Kanari’s  method 
and  the  calculations  for  the  cells  were  done  with  LiCoCL  as 
reported  by  Kanari  et  al.  [13].  A  simplified  model  was  also 
used  under  the  following  assumptions  [2,12,13]: 

1.  A  group  in  the  cell  rolls  with  both  electrodes  and  the 
separator  soaked  in  the  electrolyte  is  made  of  a 
homogeneous  material  having  effective  thermophysical 
properties. 


Table  1 

Simulation  parameters  and  their  values  for  2  Wh-class  cells 


Cell  parameters 

Values 

Diameter  (mm) 

14.2 

Height  (mm) 

50 

Initial  temperature  (K) 

298 

Environment  temperature  (K) 

298 

Heat  transfer  coefficient  (W  m-2  K~l) 

10 

250 
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transfer  coefficient  for  the  2  Wh-class  cells  was  the  mea¬ 
sured  value.  The  model  of  the  cells  was  divided  at  intervals 
of  1  mm  in  the  r  and  z  directions.  The  internal  cell  resistance 
of  each  2  Wh-class  cell  was  200  mQ,  calculated  from  the 
current-voltage  characteristics  at  the  temperature  of  298  K. 
The  initial  state  of  the  experimental  and  thermal  simulation 
in  discharge  was  0%  for  the  depth  of  discharge  maintained  at 
298  K  after  the  rated  capacity  was  charged  at  a  C/8-rate. 

4.  Results  and  discussion 

First,  we  measured  the  discharge  characteristics  of  the 
2  Wh-class  cells  in  which  LiCo02  or  LiNio.7Coo.3O2  was 
used  as  the  positive  electrode  material.  The  discharge  curves 
of  the  2  Wh-class  cells  are  shown  in  Fig.  1 .  The  cell  using 
LiNi07Co03O2  showed  a  larger  discharge  capacity  of  0.5 
than  the  0.45  Ah  displayed  by  the  cell  using  LiCo02. 

Next,  the  surface  temperatures  of  these  cells  were  mea¬ 
sured  during  discharge.  The  surface  temperatures  of  the 
2  Wh-class  cells  during  discharge  are  shown  in  Fig.  2. 
The  discharge  current  was  500  mA,  and  this  value  nearly 
corresponded  to  the  lC-rate.  Each  discharge  capacity  at  the 


Discharge  Capacity  (Ah) 


Fig.  1.  Discharge  curves  of  2  Wh-class  cells.  Negative  electrode: 
graphite-coke  hybrid  carbon;  electrolyte:  1.0  mol  P1  LiPF6-EC/DEC  (1/ 
1);  charge  current:  65  mA;  discharge  current:  500  mA. 


Depth  of  Discarge  (%) 


Fig.  2.  Surface  temperatures  of  2  Wh-class  cells  during  discharge. 
Discharge  current:  500  mA. 


258  268  278  288  298 

Temperature  (K) 

Fig.  3.  Open  circuit  voltage  vs.  temperature  of  2  Wh-class  cells  using 
LiCo02  measured  at  several  depths  of  discharge. 

lC-rate  of  the  2  Wh-class  cells  was  about  83%  of  the  charge 
capacity.  The  surface  temperature  of  the  cell  using  LiCo02 
increased  to  about  308  K  from  298  K  at  the  end  of  the 
discharge.  The  surface  temperature  of  the  cell  using  LiNi0.7_ 
Co0  3O2  increased  to  about  301  K.  The  cell  using  LiNi0.7_ 
Coo.302  showed  a  much  lower  surface  temperature  rise 
during  the  discharge  than  the  cell  using  LiCo02.  From  these 
results,  we  found  that  the  cells  with  different  positive 
electrodes  had  different  thermal  behaviors.  In  particular, 
no  significant  temperature  rise  at  the  end  region  of  the 
discharge  was  detected  for  the  cell  using  LiNi0  7Co0  302, 
but  was  seen  for  the  cell  using  LiCo02. 

Because  the  thermal  behaviors  of  the  2  Wh-class  cells 
above  were  quite  different,  we  analyzed  the  thermal  char¬ 
acteristics  of  these  cells  by  thermal  simulation.  Fig.  3  shows 
the  open  circuit  voltage  versus  temperature  of  the  2  Wh- 
class  cell  using  LiCo02  measured  at  several  depths  of 
discharge.  As  the  figure  shows,  the  open  circuit  voltage 
slightly  changes  linearly  according  to  the  temperature,  and 
the  entropy  change  “AS”  is  calculated  using  the  slope  of  the 
line.  Fig.  4  is  similar  to  Fig.  3  but  with  LiNi0.7Co0.3O2  for  the 
positive  electrode. 

From  these  results,  the  value  of  (— dE/dT)  was  calculated 
for  each  cell,  and  plotted  versus  the  depth  of  discharge  in 
Fig.  5  [14J.  The  thermal  characteristics  were  found  to 
be  quite  different  for  the  cells  using  different  positive 
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Temperature  (K) 

Fig.  4.  Open  circuit  voltage  vs.  temperature  of  2  Wh-class  cells  using 
LiNio.7Coo.3O2  measured  at  several  depths  of  discharge. 
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Fig.  5.  Relationship  between  entropy  change  and  depth  of  discharge  (a) 
LiNio.7Coo.3O2,  (b)  LiCo02. 

electrodes.  In  this  figure,  the  upper  area  of  the  line,  where  the 
value  of  (— dE/dT)  is  0.0  mV,  shows  an  exothermic  reaction 
and  the  lower  area  shows  an  endothermic  reaction  during 
discharge.  The  cell  using  LiCo02  generated  heat  all  the  time 
during  the  discharge  period.  On  the  other  hand,  the  thermal 
generation  of  the  cell  using  LiNio.7Coo.3O2  changed  accord¬ 
ing  to  the  depth  of  discharge. 

Following  this,  the  data  was  put  into  the  simulation 
equation  and  the  temperature  change  for  each  cell  during 
discharge  was  calculated.  The  results  were  compared  with 
measured  surface  temperature  values.  Fig.  6  compares  ther¬ 
mal  simulation  results  and  measured  surface  temperatures  of 
2  Wh-class  cells  during  discharge.  In  this  figure,  the  thermal 


Fig.  6.  Comparison  of  thermal  simulation  results  and  measured  surface 
temperatures  of  2  Wh-class  cells  during  discharge. 


Table  2 

Specifications  of  a  250  Wh-class  cell 


Positive  material 

LiNio.7Coo.3O2 

Negative  material 

Graphite-coke  (4/1) 

Electrolyte 

LiPF6-EC/DEC 

Diameter  (mm) 

64.0 

Length  (mm) 

296 

Weight  (kg) 

2.2 

Table  3 

Simulation  parameters  and  their  values  for  a  250  Wh-class  cell 


Cell  parameters  Values 


Diameter  (mm)  64.0 

Height  (mm)  296 

Initial  temperature  (K)  298 

Environment  temperature  (K)  298 

Heat  transfer  coefficient  ( W  m  2  K  1 )  7 


simulation  results  show  a  good  agreement  with  the  measured 
temperatures.  In  addition,  it  is  clear  that  the  cell  using 
LiNio.7Coo.3O2  has  a  smaller  heat  generation  than  that  using 
LiCo02. 

Next,  the  thermal  simulation  model  was  applied  to 
250  Wh-class  large  cells.  A  250  Wh-class  cell  is  shown  in 
Fig.  7  and  the  specifications  of  the  250  Wh-class  cell  are 
shown  in  Table  2.  The  cell  parameters  and  their  values  are 
listed  in  Table  3.  The  heat  transfer  coefficient  for  the 
250  Wh-class  cells  was  the  measured  value.  Like  above, 
the  model  of  the  cells  was  divided  at  intervals  of  1  mm  in  the 
r  and  z  directions.  The  internal  cell  resistance  of  a  250  Wh- 
class  cell  was  8  mQ,  calculated  from  the  current-voltage 
characteristics  at  the  temperature  of  298  K.  LiNio.7Coo.3O2 
was  used  for  the  positive  electrode  material,  and  the  other 
materials  were  the  same  as  those  used  for  the  2  Wh-class 
cells  except  for  an  aluminum  vessel.  The  change  in  voltage 
and  the  surface  temperature  were  measured  during  the  initial 
charge-discharge.  The  charge  current  was  8.8  A  (which 
nearly  corresponded  to  the  C/8-rate)  and  the  discharge 
current  was  70  A  (which  nearly  corresponded  to  the  1C- 
rate). 

Fig.  8  shows  the  charge/discharge  curves  of  a  250  Wh- 
class  cell.  The  discharge  capacity  of  this  cell  was  about  97% 
of  the  charge  capacity.  Fig.  9  shows  the  surface  temperature 


Fig.  7.  250  Wh-class  cell. 
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Fig.  9.  Thermal  simulation  results  and  measured  surface  temperatures  of  a 
250  Wh-class  cell  during  discharge.  Charge:  C/8  (8.8  A);  discharge:  1C 
(70  A). 


of  a  250  Wh-class  cell  measured  during  discharge  and  the 
thermal  simulation  results.  The  thermal  simulation  results 
closely  match  the  actual  measurements.  But  the  measured 
surface  temperature  values  in  the  latter  half  of  the  measure¬ 
ment  period  are  lower  than  the  thermal  simulation  values. 
We  believe  that  this  difference  resulted  from  the  decrease  in 
the  resistance  of  the  cell  by  the  heat  generation.  A  future  task 
will  therefore,  be  to  incorporate  the  dependence  of  the 
resistance  on  the  temperature  and  depth  of  discharge  into 
the  thermal  simulation. 


5.  Conclusions 

The  results  in  this  paper  are  summarized  as  follows: 

1.  The  thermal  simulation  results  of  a  2  Wh-class  cell 
showed  good  agreement  with  measured  temperature 
values. 


2.  The  heat  generation  of  a  cell  using  LiNio.7Coo.3O2  was 
less  than  that  using  LiCo02  during  discharge.  This 
difference  was  caused  by  the  difference  in  the  change  of 
entropy.  LiNio.7Coo.3O2  was  found  to  be  a  suitable 
positive  electrode  material  for  large-scale  cells  because 
of  its  small  heat  generation. 

3.  A  250  Wh-class  cell  using  LiNio.7Coo.3O2  was  con¬ 
structed,  and  its  simulation  results  were  in  accord  with 
measured  temperature  values.  The  simulation  was  found 
to  be  efficient  in  estimating  the  thermal  characteristics 
of  both  small-sized  cells  and  large-scale  cells. 
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